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Abstract: We evaluate the probability of future neutrino oscillation facilities to 
discover leptonic CP violation and/or measure the neutrino mass hierarchy. We study 
how this probabihty is affected by positive or negative hints for these observables to 
be found at T2K and NOuA. We consider the following facilities: LBNE; T2HK; and 
the 10 GeV Neutrino Factory (NFIO), and show how their discovery probabilities 
change with the running time of T2K and NOz/A conditioned to their results. We 
find that, if after 15 years T2K and NOz/A have not observed a 90% CL hint of 
CP violation, then LBNE and T2HK have less than a 10% chance of achieving a 
5cr discovery, whereas NFIO still has a ~ 40% chance to do so. Conversely, if T2K 
and NOi/A have an early 90% CL hint in 5 years from now, T2HK has a rather 
large chance to achieve a 5a CP violation discovery (75% or 55%, depending on 
whether the mass hierarchy is known or not) . This is to be compared with the 90% 
(30%) probability that NFIO (LBNE) would have to observe the same signal at 5a. 
A hierarchy measurement at 5a is achievable at both LBNE and NFIO with more 
than 90% probability, irrespectively of the outcome of T2K and NOz/A. We also find 
that if LBNE or a similar very long baseline super-beam is the only next generation 
facility to be built, then it is very useful to continue running T2K and NOz/A (or 
at least T2K) beyond their original schedule in order to increase the CP violation 
discovery chances, given their complementarity. 
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1. Introduction 

A coherent picture of leptonic mixing is gradually emerging from the experimental 
data, much as in the case of quark mixing: the three Standard Model lepton gener- 
ations mix via the unitary Pontecorvo-Maki-Nakagawa-Sakata [1-5] matrix ?7pmnS) 
the analogue of the CKM matrix Vckm that governs hadronic mixing. This matrix 
can be parametrized in terms of three rotation angles, 9i2, Ois and 623, that have 
been first measured using solar [6-13], reactor [14-16] and atmospheric [17] neutrino 
data, respectively. Their values have been further constrained using long-baseline 
reactors [18] and accelerator-based neutrino beams [19-24]. The most recent global 
fits of neutrino data give the following best-fits and la errors for the leptonic mixing 
angles and the neutrino squared mass differences [25] (see also Refs. [26,27]): 

r sinHu = 0.3021^^2 f Amii/10-5 = 7. 50 trolley 

sin^ ^13 = 0.0227l°:°°23 J AmlJlQ-^ = 2.473l°:°™eV' (NH) 



[ AmlJlO-' = -2.427+%°6t eV^ (IH) 



sin^^,3 = 0.413 70.5941^22 

where the twofold ambiguities related to allowed values of the ^^23-octant and of the 
neutrino mass hierarchy are clearly shown. These values strongly differ from the 
corresponding values obtained for the mixing angles in the CKM matrix [28], when 
the latter is also parametrized in a similar way. It is well known that the CKM matrix 
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can be approximately parametrized as a small perturbation of the identity matrix, 
while the PMNS matrix rather displays entries which are all of similar order. The 
different patterns of quark and lepton mixing add to the puzzling flavor problem, i.e. 
the apparent arbitrariness of the Yukawa couplings of the Higgs boson with fermions 
belonging to different generations, that eventually induce fermion masses with large 
hierarchies after the spontaneous breaking of the electroweak symmetry. The second 
important difference between the CKM matrix and the PMNS matrix is that in the 
former case we know that the matrix is complex, i.e. that a non- vanishing CP- 
violating phase is needed in order to accurately reproduce the hadronic flavor mixing 
data. The quark CP-violating phase is 5ckm ~ 7 = (66 ± 12)° (see, e.g., Refs. [29] 
and [30]). A straightforward question is, therefore, if the leptonic mixing matrix is a 
complex matrix, as well as Vckm- 

The task of establishing whether (5pmns is CP-violating (and possibly measuring 
it with a precision comparable to that of 5ckm) is the main goal of the present 
generation of long baseline neutrino experiments. T2K [31] is presently running, 
and the NOi/A experiment [32] is under construction and is expected to start data 
taking soon (in May 2013). Even if these experiments will be able to probe CP 
violation, their combined significance will remain below the 3cr evidence unless they 
run for an extended period of time, see for instance Refs. [33,34]. Therefore, a next 
generation of long baseline neutrino experiments will most likely be needed if we 
want to establish leptonic CP violation at the 5a level. A second task that these 
experiments should address is the determination of the neutrino mass spectrum, i.e. 
if Am|]^ > (normal hierarchy) or Ato|^ < (inverted hierarchy). However, due to 
the rather large value of ^13 recently discovered, the measurement of the neutrino 
mass hierarchy may be accessible to several alternative strategies to a long baseline 
neutrino experiment, such as atmospheric neutrino experiments [35-39] or a reactor 
experiment with a medium baseline [40-47]. 

In this paper we compare the proposed new facilities to be built in the following 
ten to twenty years based on their chances to establish a CP-violating phase 5pmns 
(from now on S, for simplicity) at a given CL and, secondarily, on their chances to 
pin down the correct neutrino mass hierarchy. For this purpose, we compute the 
probability that a given new facility has to discover a CP-violating 5 (or a given 
mass ordering) at 5(7, conditioned to the presence or absence of a hint (at the 90% 
CL for CP violation and at the 3a level for the mass hierarchy)^ for the combination 
of T2K+N0z/A. We define a condensed performance indicator, the gain fraction G 
(and G), corresponding to the conditional probability that a next generation facility 
would reach the discovery limit over an observable, given that T2K-|-N0i/A have (or 

'^Whilc higher CL can prove chahenging to the statisticaUy-hmitcd T2K+N0z^A, they will be 
able to provide 90% hints over a large portion of the parameter space. Therefore, while we will 
require that the new experiment reaches the discovery limit, for T2K-|-N0i/A only a hint at the 
90% CL will be required. 
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have not) previously observed a hint of it. The gain fraction thus represents the 
"gain" in performance that the new facihty can provide given the prior results from 
T2K and NOz/A. This approach allows for a quantifiable betting strategy between 
the proposed new facilities, particularly important at this time of decision-making 
in neutrino oscillation experiments. 

The paper is organized as follows. In Sect. 2 we introduce the statistical approach 
and the terminology used in the rest of the paper before describing the technical 
details of the included facilities in Sect. 3. We then present our results in Sect. 4, 
before we summarize and give our conclusions in Sect. 5. 

2. The statistical method 

We address the question of how much discovery potential can be gained by building a 
new facility given the prior information available at the time of its running from T2K 
and NOz/A. If a significant hint for CP violation or a given mass hierarchy is already 
present, then a modest upgrade to increase the statistics should be able to provide 
the desired CL. Conversely, if no hint is found, a more demanding facility is needed 
to significantly increase our discovery potential. In other words, depending on the 
outcome of T2K and NOi^A, can we quantify the additional gain a next generation 
facility would provide? 

We will make use of some simple statistical tools. First of all, we compute the 
probability that a given facility combined with T2K+N0z/A reaches a desired CL 
(outcome A) for a given observable (cither discovery of CP violation or of the correct 
mass hierarchy) as a function of the true value of 6, P{A\6). In order to compute 
this probability we take the following steps: 

1. First, the expected event rates (including signal and backgrounds) are com- 
puted for each energy bin and oscillation channel of the facility under study, 
for an input value for S; 

2. A Monte Carlo is then used to generate 1000 "realizations" of the experiment, 
by randomly generating numbers for the different energy bins and channels ac- 
cording to a Poisson distribution with mean value equal to the one computed 
in step 1. In order to do this, a new function has been added to the Monte- 
CUBES free software [48] that produces a random distribution of events at a 
given experiment, provided its expected event rates; 

3. For each realization of the experiment, we evaluate the CL with which CP 
violation or the correct mass hierarchy is favored. Then, for a given value of 5, 
we count the fraction of times the target CL was reached, that is, the fraction 
of times for which A is fulfilled, thus estimating the probability P{A\5). 
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The probability P{A\6) is the starting point of our statistical analysis. In particular, 
instead of computing the expected achievable CL for a given observable {e.g. the 
capability to observe a CP-violating signal) at a given facility, we rather fix a target 
CL (typically 90%, 3a or 5a) and ask what is the probability with which the facility 
will reach the target. This procedure is slightly different from what is commonly 
done in the literature, where the most widely used performance indicator is rather 
the discovery potential for a given observable, defined as the CL that a facility is 
expected to provide for that observable as a function of S. In this case, the expected 
number of events for the facility is computed for each value of S and the CL with 
which CP conservation or the wrong mass hierarchy would be disfavored for that 
value of 6 is evaluated. It has been shown (see Ref. [49]) that the probability P{A\6) 
of the facility actually reaching that CL is around 50%. Indeed, upward or downward 
random fiuctuations will increase or decrease the final CL in a particular realization 
of the experiment. 

To address the question of how much we can improve on the results from T2K 
and NOi/A we will also make use of the joint probability P{A,B\S), representing 
the probabihty that the new facihty combined with T2K-|-N0i/A reaches a given CL 
(outcome A) and that the combination of T2K+N0i^A alone also observed a signal at 
some CL for the same observable (outcome B)^, for a given S. It should be noted that 
the CL required for outcomes A and B do not need to be the same. The computation 
of the joint probability is done following similar steps to 1-3. We first compute the 
expected event rates for both T2K and NOz/A and the next generation facility and 
generate 1000 different reafizations of them. Then, for each reafization, we record the 
CL with which T2K and NOi^A, on their own as well as in combination with the new 
facility, are able to disfavor CP conservation or the wrong mass hierarchy. Finally, 
we count how many times T2K and NOz/A were (un)able to reach the desired CL 
while the combination with the new facility reached its desired CL. This procedure 
is repeated for all values of S. 

We are interested, in particular, in the joint probability that T2K+N0z/A find 
(or not) a hint at the 90% CL for a given value of 6 {B or B, respectively) and 
that their combination with a new facility reaches 5a for the same 6 and observable 
(A). In the case of the mass hierarchy we will raise the CL of B to the 3a level, 
since T2K-FN0i/A will most hkely have a 90% hint. The joint probability P{A, B\d) 
defined this way is a measure of the complementarity between T2K-|-N0i/A and the 
new facihty for that particular value of S. Indeed, if P{A, B\S) is high it means that 
there is a high probability that T2K-|-N0i'A will not find a hint at that value of S 
while the probability of making a discovery of their combination with the new facility 
is also high. Clearly, for facilities whose best chances to measure S take place for the 

^We will also let B represent the complement oi B, i.e., that T2K-|-N0i/A do not reach a given 
CL. We then have P{A, B\6) + P{A, B\6) = P{A\5). 
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Event 


Definition 


A 


The combination of T2K+N0z/A witli a future facility provides a 
discovery of CP violation or the mass hierarchy at 5o". 


B 


T2K+N0z/A provide a signal for CP violation (or mass hierarchy) 
at 90% CL (3(7). 


B 


T2K+N0z/A do not provide a hint. The complement of B. 



Table 1: The definitions of the events A, B, and B. Cited for hack reference. 



same region of parameter space as for the combination of T2K and NOi/A, the joint 
probability P(A, B\S) will be small, since for favorable values of 5 a positive result at 

the 90% CL would have been likely at T2K+N0z/A. On the other hand, for a facility 
whose best sensitivity to 5 is not totally overlapping with that of T2K+N0z/A, we 
expect a rather high joint probability in the regions where they do not overlap. In 
order to simplify for the reader to recollect our definitions of A and B, these can be 
found in Tab. 1. 

A very interesting (and condensed) quantity to parametrize the margin for im- 
provement over the results of T2K-|-N0i/A by a new facility is the conditional prob- 
ability P{A\B) — P{A, B)/P{B). This is the probabihty of achieving a 5(7 discovery 
at the new facility provided a 90% CL or 3(T signal has been previously obtained at 
T2K+N0z/A. For brevity, we will call this quantity the gain fraction G. Similarly, we 
can define the gain fraction assuming the opposite outcome: G = P{A\B), i.e., the 
probability of achieving a 5cr discovery at the new facility provided no 90% CL or 3a 
signal has been previously seen at T2KH-N0i/A. Notice that we have defined the gain 
fractions in a way that does not depend on S. Indeed, the outcome from T2K-|-N0i/A 
will provide a probability distribution for the possible values of 5, P{S\B), that can 
be used as prior information for the measurements of the new facility. Integration 
can be performed over 5 as follows^: 

(...) 

We can see that there is a very close relation between the joint probability P{A, B\S) 
and the gain fractions. Thus, the joint probability is not only useful as a measure 
of the complementarity of a new facility with respect to T2K+N0z/A, but can also 
help to understand the behaviour of the gain fraction. For this reason, example plots 
of the joint probabihty for the CP and mass hierarchy discovery as a function of S 
will also be shown in Sect. 4. We will also show the dependence of the gain fractions 
G and G on the running time of T2K-|-N0i/A. Prom the above equation, it is clear 



^Using that P{S\B) = P{B\5) x P{5)/P{B) according to Bayes' theorem. We also assume a flat 
distribution for 5, P{5) = l/(27r), prior to the measurement from T2K+N0z/A. 
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2.5 


Wu 


22.5 
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0.2-0.7 




Var. 


NOi/A 


2 


810 


0.8° 


TASD 


3-14 


0-0.7 




Var. 


NFIO 


6 


2000 




MIND 


100 




7x1020 


(10,10) 


LBNE 


3.0 


1290 




LAr 


10-33 


0.8 




(5,5) 


T2HK 


0.6 


295 


2.5° 


WC 


560 


1.66 




(1.5,3.5) 



Table 2: Main features of the setups considered in this work. From left to right, the 
columns list the energy at which the number of neutrino events peak (in GeV), the length of 
the baseline (in km), the off-axis angle in degrees, the detector technology (Magnetized Iron 
Neutrino Detector, Liquid Argon, Water Cerenkov or Totally Active Scintillator Detector), 
the fiducial mass of the detector, the beam power or the number of useful decays per year, 
and the number of years that the experiment is assumed to run per polarity. The label 
"Var. " indicates that the running times for T2K and NOvA are variable, see text for 
details. 

that the gain fraction can also be obtained in a graphical way: it is simply the area 
under the curve for P{A,B\5) divided by the area under P{B\5), thus representing 
the fraction of parameter space remaining after T2K-|-N0i/A that is probed by the 
new facility. For zero running time of T2K-|-N0i/A the gain fraction is, therefore, 
intimately related to the CP-fraction that is usually discussed in the literature. 

Finally, it should be noted that B and B are oversimplifications of the results 
from T2K+N0z^A. Indeed, the outcome from T2K+N0z/A will not be limited to the 
presence or absence of a hint, but will rather be some measured event distributions 
that would translate in a given probability distribution for 5 when trying to fit 
different values to data, P(5|data). It is this probabihty distribution which should be 
used as condition B for the prior when computing the gain fraction. However, lacking 
the results from T2K-|-N0i/A we decided to adopt two representative scenarios to 
present our results: the presence or absence of a 90% CL or 3(7 signal. 

3. Setups 

The performance of new facilities will be compared to (and combined with) those 
attainable from the combination of T2K and NOz/A, assuming variable running times. 
The main features for these two experiments are summarized in Tab. 2. Notice that 
the beam power of T2K has not reached its goal yet and it is not expected to do so 
in the next couple of years. In the case of NOi/A, both the fiducial volume of the 
detector as well as the beam power are expected to increase progressively within the 
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Figure 1: Integrated exposure as a function of time for NOvA (left panel) and T2K (right 
panel). For NOvA, the running time is split equally between neutrino and antineutrino 
modes. The shown integrated exposure corresponds to the total between the two polarities. 
For T2K, the exposure for each polarity is assumed to change differently with time, as seen 
in the right hand panel. The integrated exposure in this case corresponds to each polarity 
separately. The vertical lines correspond to year 2020. 



first few montlis of data taking of the experiment, which is expected to start in May 
2013 [50]. Several assumptions have been made in this work in order to perform as 
reahstic as possible an estimate of the statistics achievable at each experiment for a 
given date. Fig. 1 shows our assumptions for the integrated exposure for T2K and 
NOi^A as a function of time. The integrated NOz/A exposure has been extracted 
from Ref. [50]. For T2K we have followed Ref. [51], assuming that the experiment 
starts in 2013 at 200 kW, goes up to 300 kW in 2014, and finally reaches 700 kW 
in 2017. In addition, we have assumed that during the first 5 years it is only run in 
neutrino mode, but that after that date the running is equally split between the two 
polarities. A vertical line corresponding to year 2020 has been added to Fig. 1 for 
reference. 

The simulation of the T2K and NOz/A experiments is largely based on Ref. [52]. 
The T2K simulation is exactly the same, except for the integrated exposure which 
has been varied according to Fig. 1. The situation is a bit different for NOz/A. In 
view of the large value of ^13 the NOz/A collaboration have recently modified their 
analysis techniques [50]. In particular, some of the cuts have been relaxed so that 
larger efficiencies are now expected, at the price of larger backgrounds. Therefore, we 
have used the same signal and background rejection efficiencies as in Ref. [53], which 
are in agreement with those in Ref. [50]. Finally, we are using matrices to migrate 
NC backgrounds to low energies, unlike in Ref. [52] where migration matrices were 
not available for this purpose. 

Alongside with T2K and NOz/A, we consider in this paper a few representative 
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facilities to analyse under the new performance indicators described in Sect. 2. We 
will concentrate in particular on three of the most discussed setups for the next 
generation of neutrino oscillation facilities: two super-beams, T2HK and LBNE, 
and the "low energy" [54-58] Neutrino Factory (NFIO), which is considered as the 
ultimate neutrino oscillation facility but also a longer term project. 

The technical details of the simulations for these three facilities are summarized 
in Tab. 2: 

1. T2HK 

This setup uses the same baseline (L = 295 km) and accelerator as T2K (al- 
beit, with a higher power of 1.66 MW) aiming at a new detector, Hyper- 
KamiokaNDE (HK), with a fiducial volume of 560 kton of water [59]. To 
compensate for the different vN and VN CC cross-sections, this setup is fore- 
seen to run for 1.5 years with tt"*" and for 3.5 years with tt". The collaboration 
is also exploring the possibility to run the experiment for 10 years (3 years with 
TT"*" and 7 years for vr^) but with a beam power of 750 kW (see, for instance, 
Refs. [60,61]). We have checked that the results using this configuration are 
very similar, and therefore will only show those corresponding to the configu- 
ration described in the HK letter of intent, Ref. [59]. 

2. LBNE 

The proposal consists of a new accelerator to be built at FermiLab with a 
design power of 0.7 MW, aiming at a distant, on-axis, 33 kton Liquid Argon 
(LAr) detector (in a first phase, only 10 kton of LAr are considered) located 
at L = 1290 km from the source"^. Five years of data taking for tt"'" and tt" 
are considered. The simulation details have been implemented according to 
Ref. [66]. 

3. NFIO 

After the measurement of a large ^^13, a rc-cvaluation of the Neutrino Factory 
design [67] has identified the optimal parent muon energy for this facility at 10 
GeV (see also Ref. [68]). The new design includes a 100 kton magnetized iron 
neutrino detector (MIND, [69]), located at approximately L — 2000 km from 
the source. Both /x+ and would run at the same time inside the decay ring 
for a total of 10 years. The detector simulation details have been implemented 
according to Ref. [69] . Migration matrices for both signal and backgrounds have 
been used [70] , and the so-called r-contamination [71-74] has been included in 
the analysis for both appearance and disappearance signals. 

^European alternatives to these two facilities with somewhat similar physies performances have 
also been proposed. The SPL [62, 63] or ESS proposals [64] present similar characteristics and 
performance to T2HK, while LBNO [65], with a smaller LAr detector, longer baseline and higher 
energy, has a similar performance to LBNE. 
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Systematic errors are expected to have a large impact on the resuhs for long 
baseline experiments, given the large value of 6i^. In Ref. [75] a detailed analysis 
of systematics and their impact on long baseline oscillation experiments was per- 
formed. It was shown that the effect on LBNE as well as the NFIO setups is small^ 
if a fully correlated analysis is performed. In the case of T2HK, however, the re- 
sults present large variations depending on the particular values for the systematic 
uncertainties and, more importantly, on the particular implementation used. In the 
present work, systematic uncertainties have been implemented as overall normaliza- 
tion errors, which are uncorrelated between signal and backgrounds but correlated 
among all backgrounds contributing to the same channel. Even if this does not cor- 
respond to the real physics case, it was shown in Ref. [75] that an appropriate choice 
of the normalization errors would yield similar results to those obtained with a more 
sophisticated systematics implementation. Thus, we have chosen the following nor- 
malization errors for the signal and backgrounds in each experiment: 1.5% and 10% 
for the NFIO, 3% in both cases for LBNE, and 5% and 10% for T2HK and NOuA. 
For T2K, on the other hand, 2% and 5% normalization errors have been used for 
the signal and background event rates, plus a calibration (tilt) error of 1% and 5%, 
respectively, following Ref. [52]. 



4. Results 

In this section we apply the statistical method described in Sect. 2 to compare the 
different facilities defined in Sect. 3, under the hypotheses that the combination of 
the data taken by T2K and NOi/A, after a given number of years of running time, 
either does or does not show a hint of a CP-violating phase S (or of a given hierarchy) 
at the 90% CL (3(7). We restrict to this CL due to the statistical limitations of both 
experiments. 

4.1 CP violation 

Fig. 2 shows the results for T2K-|-N0i/A assuming that they will take data until 
2020 {i.e., close to the initially foreseen running time). The probability of finding 
CP violation at 90% CL for the combination of T2K-|-N0z^A is represented in Fig. 2 
(left) by the gray-shaded area. In the same panel, the lines show the probability of 
finding CP violation {i.e., to distinguish a CP-violating phase 6 from and tt) at 5a 
for the different new facilities described in Sect. 3, in combination with T2K+N0z/A 
data. Solid black lines correspond to NFIO; dotted blue lines to T2HK; dot-dashed 
red lines to LBNE- 33 kton; and dashed green lines to LBNE- 10 kton. As explained in 

^Note that no uncertainties on the shape of the neutrino flux or cross section were considered 
in Ref. [75], though. If this is not the case the situation may be very different, see for instance 
Refs. [76-78]. 
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Sect. 2, this probability is obtained by simulating different realizations of a given ex- 
periment 1000 times for each input value of 6 and checking whether in that particular 
realization the facility is able to exclude CP-conserving values of S at the desired CL. 
For T2K+N0z/A and for T2HK, the probability of finding CP violation is strongly 
asymmetric, reaching higher probabilities for negative 6. This is a consequence of 
the sign degeneracy that leads to allowed CP-conserving regions for positive S when 
measuring — )■ z/g oscillations [79-81]. 
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Figure 2: Left panel: probability of finding CP violation (CPV) at the 90% CL for the 
combination T2K+N0iyA (gray-shaded area); probability of finding CP violation at 5a for 
the new facilities in combination with the results obtained at T2K+N0uA (lines). Right 
panel: joint probability that a new facility in combination with T2K+N0vA observes CP 
violation at 5a and that T2K+N0vA alone have not seen a hint of it at 90% CL for the 
same value of 5. In all cases, it is assumed that T2K+N0uA will stop data taking on 2020. 



To explain our results let us consider now, for example, the black solid line in 
Fig. 2 (left), that represents the probability of finding CP violation at So" (condition 
A) for^ the NFIO, P{A\5). Two flat regions can be seen for 6 G [-150°, -30°] and 
[30°, 150°], for which the probability to measure a CP- violating phase at 5cr at the 
NFIO is P{A\5) ~ 1. Depending on 5, it is more or less likely that T2K-|-N0z/A 
will have already found a 90% CL hint by themselves. For maximally CP-violating 
values of 5, in particular, there is a fair chance that they will have found a hint 
by 2020 and, thus, condition B will not be satisfied. This is reflected in that the 
black solid line in Fig. 2 (right), which represents the joint probability P{A, B\6), 
is reduced for \6\ = 90°. For S = —90°, where T2K-|-N0z/A have a peak in the 
probability of finding CP violation at the 90% CL, we see that P{A\ - 90°) ~ 1 
whereas P{A,B \ —90°) ~ 0.3. On the borders of the region of maximal probability of 
finding CP violation for T2K-|-N0i/A we see high spikes for the NFIO joint probability. 

^We have checked that the combination of the NFIO data with T2K+N0i'A is irrelevant, due 
to the low statistics of the latter with respect to the NFIO. 
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These appear because the flat region that we observe in Fig. 2 (left) for NFIO is much 
broader than the region where the probabihty of flnding CP violation for T2K+N0z/A 
(the gray-shaded area) is maximal. This shows that, in those spikes, NFIO can probe 
for CP violation with very high significance, while T2K-|-N0i^A have a hard time 
finding a hint. Thus, the gain factor G for NFIO will receive significant contributions 
from those areas of the parameter space. 




-150-100 -50 50 100 150 



<5[°] 

Figure 3: Joint probability that NFIO observes CP violation at So" for a fixed value of 6 
and that TBK+NOvA have not been able to see a hint of it at the 90% CL for the same 
value of 5. Results are shown assuming that T2K+N0vA stop data taking in 2015, 2020, 
2025 or 2030, as indicated in the legend. 

It is interesting to observe how the joint probability P{A, B\S) is affected as the 
running time of T2K+N0z/A is increased. This is illustrated in Fig. 3, where we 
show the joint probability that NFIO discovers CP violation at 5(j (condition A) and 
T2K-|-N0z/A have not seen a hint of it at the 90% CL if they stop data taking on 2015, 
2020, 2025 or 2030 (condition B). It can be clearly seen that the probability that 
both conditions are fulfilled for maximal CP-violating values of 5 (5 ~ ±90°) rapidly 
decreases with increasing T2K-|-N0i/A exposures, whereas spikes at 5 ~ ±30°, ±150° 
are only mildly affected (showing that, for NFIO, some areas of the parameter space 
remain promising even after a negative result from T2K+N0z/A). The reason for 
this is that, as the data taking period increases for T2K+N0z/A, it is more and 
more unlikely that they do not observe a hint for CP violation at the 90% CL if S is 
maximally CP-violating. 

The same effect is illustrated in Fig. 4 for all facilities, which shows the same 
probabilities as in Fig. 2 but assuming that T2K+N0i/A run for 10 additional years, 
until 2030. As it can be seen, the probabihty of not finding a positive 90% CL hint 
for S ~ ±90° if T2K+N0i/A keep running until 2030 is very small. This in turns 
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Figure 4: Same as in Fig. 2 assuming that T2K+N0uA will run until 2030. 



implies that, provided a negative result from T2K+N0i/A, CP-violating values for 
5 are much less likely and the new facility will have smaller chances of finding CP 
violation. The only facility that retains a non-negligible joint probability after this 
exposure is NFIO. This fact is better captured by the gain fractions G and G defined 
in Sect. 2, which correspond to the conditional probability that the new facility will 
detect CP violation at 5cr, provided that T2K-fN0z/A have previously seen (G) or 
have not seen [G) a hint at the 90% CL. 

We show in Fig. 5 bands for G (left panel) and G (right panel) as a function of 
the year when T2K-|-N0z/A stop data taking. The lower and upper limits of each 
band correspond to two possible scenarios: (1) the neutrino mass hierarchy is not 
known at the time when the data of T2K-|-N0i/A are analyzed, and (2) the neutrino 
mass hierarchy is known at Scr at that time, assuming that other experiments have 
been able to measure it. The gain fractions obtained in case 1 are depicted by thick 
lines, whereas those corresponding to case 2 are shown as just the end of each band. 
As shown in the figure, G always increases when the additional information about the 
mass hierarchy is added, while the situation is the opposite for G (with the exception 
of T2HK). In case the mass hierarchy is measured during a certain year, it suffices 
to consider the band edge corresponding to case 2 to take into account the impact 
of this information on the gain fractions. It can be seen that previous knowledge on 
the hierarchy plays a very important role for the performance of T2HK, while it is 
less significant for other new facilities (albeit, G is generically more affected than G). 
Color codes for each facility are shown in the legend. 

Consider first Fig. 5 (left), which shows G as a function of the year when 
T2K-fN0z/A stop data taking. On general grounds, if a positive hint of a CP- 
violating phase is found at the 90% CL at T2K+N0z/A the gain fraction will increase, 
since CP-violating values of 5 become favored by data. When increasing the expo- 
sure, if we keep requiring only a 90% CL hint for CP violation, the gain fraction will 
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Year Year 



Figure 5: The gain fractions G (left panel) and G (right panel) for CP violation defined 
as in Sect. 2, i.e. the probability that the new facility will detect a CP-violating signal at 
5a provided that T2K+N0uA have seen or not a hint of it at 90%CL irrespectively of the 
value of 5. In both cases results are shown as a function of the year when T2K+N0vA stop 
taking data. Bands show how the gain fractions change if we assume no previous knowledge 
of the mass hierarchy (thick edges) or that the mass hierarchy is known by the time the 
new facility analyses its data (thin edges). 



tend to a plateau. If we were to use real T2K+N0z/A data instead of just requiring a 
90% CL hint, the CL at which CP violation would be favored should rather increase 
with the exposure (assuming that 6 is indeed CP-violating) and hence the gain frac- 
tion would continue rising instead of reaching a plateau. In Fig. 5 (left), when the 
mass hierarchy is known, the plateau is clearly visible for NFIO and T2HK. However, 
when the data from T2K+N0z/A are complementary to that of the new facility and 
carry significant statistical weight, the gain fraction will continue to rise with the 
exposure of T2K-fN0i/A, since they play an important role in reaching the 5a goal. 
This behaviour is clearly seen for both LBNE setups, which show a monotonously 
increasing G with T2K+N0i/A running time, reaching ~ 10%(~ 40%) by 2030. For 
these facilities, the combination with more and more T2K+N0i^A data will increase 
the overall probability to measure CP violation at 5cr. 

Notice also that the behavior of G for NFIO, LBNE-10 kton and LBNE-33 kton 
is practically unaffected by the previous knowledge of the neutrino mass hierarchy 
(represented by upper edges on the bands), and when the mass hierarchy is known the 
gain fraction G only slightly increases. This is due to the fact that, at T2K-|-N0z/A, 
positive (negative) CP-violating values of 6 with a normal (inverted) hierarchy are 
very degenerate with CP-conserving values of 6 and an inverted (normal) hierarchy. 
It can then happen that, for CP-conserving values of 6, the CP- violating degenerate 
solution is actually favored at the 90% CL from T2K+N0z/A, thus leading to a fake 
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hint. These cases will not lead to a discovery at any of these new facilities that 
have very good sensitivity to the mass hierarchy and will therefore be able to rule 
out these fake hints. Conversely, if the mass hierarchy was known, these fake hints 
would not take place at T2K+N0i/A and a positive hint will more likely imply a 
truly CP- violating value for S, thus increasing the chance of discovery at the new 
facility, that is, the gain fraction. 

A completely different behavior is shown by T2HK (blue band): in the case 
the mass hierarchy is not known we see that G decreases with time, reaching a 
plateau at G ~ 0.45 by 2022. This behavior can be understood as follows: if a CP- 
violating 6 is detected for a very low exposure at T2K+N0z/A, then 6 is probably 
very near to 6 = —90° for a normal hierarchy or 6 = 90° for an inverted hierarchy, 
a value for which the combination of T2K-|-N0i/A data with those of T2HK will 
have a very high chance to measure CP violation at 5a. On the other hand, if CP 
violation is detected at 90% CL for increasing exposure, different values of S can 
give rise to the 90% CL hint. For many of these values, the sign degeneracies will 
be locahzed at CP-conserving values. And, while a high exposure from T2K-|-N0i/A 
can disfavor the CP-conserving degeneracies at the 90% CL wc arc demanding, the 
combination of T2K+N0z/A and T2HK cannot do the same at 5a and hence cannot 
discover CP violation. This explains the decreasing gain fraction of T2HK when 
increasing T2K-t-N0z/A exposure if the mass hierarchy is not known. This problem 
is of course solved if the neutrino mass hierarchy were known when T2K-|-N0i/A 
data are analized, as shown by the upper line of the T2HK blue band, for which 
we recover the expected increasing G with time. Notice that also T2HK reaches a 
plateau by 2020 (as the NFIO), with G ~ 0.7. Running T2K+N0i/A beyond this 
date, therefore, does not increase the discovery chances, unlike the case of the LBNE 
setups, which keep growing after 2020, albeit more slowly. 

Let us consider now Fig. 5 (right), showing (5 as a function of the T2K+N0z/A 
total running time. On general grounds, when no hint for CP violation is found at 
T2K-|-N0z/A as the exposure is increased, it is more and more likely that 6 is close to a 
CP-conserving value, thus making a discovery of CP violation impossible. Therefore, 
we expect that G will monotonously decrease with T2K-|-N0i/A total running time. 
Notice that the knowledge of the neutrino mass hierarchy will have a negative impact 
on the G gain fractions (with the exception of T2HK), since CP- violating regions 
disfavored at the 90% CL by T2K+N0z/A with respect to CP-conserving ones with 
opposite hierarchy could actually turn out to be the true solution and lead to a 5a 
discovery at the new facility. These fake CP-conserving solutions cannot take place 
when the mass hierarchy is known and hence a negative hint from T2K-|-N0z/A 
is more likely to imply a CP-conserving 6 in this scenario, decreasing the value of 
G. In the case of T2HK, on the other hand, previous knowledge of the hierarchy 
would increase G, since this facility is unable to measure the hierarchy by its own at 
the required 5a level and, therefore, its performance is significantly affected by sign 
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degeneracies. When we study G for T2HK, then, we notice that in case 1, represented 
by the lower edge of the blue band, we get G ~ 0.2 (0.05) for T2K+N0z/A running 
until 2015 (2030). On the other hand in case 2, represented as for G by the upper 
edge of the band, G goes from G ~ 0.45 — )■ 0.07 in the same time interval. NFIO and 
LBNE-33 kton show the expected decreasing G with T2K+N0i/A running time: for 
NFIO we have G ~ 0.7 — t- 0.4 in case 1 and G ~ 0.7 — )■ 0.35 in case 2, whereas for 
LBNE-33 kton we have G ~ 0.13 ^ 0.10 in case 1 and G ~ 0.13 in case 2. The 
case of LBNE-10 kton is a bit different: the statistical power of this facility is limited 
to its capability of measuring the neutrino mass hierarchy and, therefore, G increases 
with time since the new facility increases the ability of T2K+N0z/A actually making 
a discovery (by ruling out fake CP-conserving sign degeneracies). On the other hand, 
if previous knowledge of the hierarchy is at hand, the facility has no capability to 
discover CP violation by its own. 

We can see that if T2K+N0i/A see no hint of leptonic CP violation with in- 
creasing exposure, a clear gap between NFIO and the rest of facilities opens: we can 
fairly say that the super-beam-based experiments have a very low or null chance to 
measure CP violation, whereas NFIO would still have a 40% chance to measure it. 
On the other hand, if T2K-|-N0t^A have a positive signal already, T2HK appears 
to be an interesting alternative to the neutrino factory, in particular if some other 
experiment is able to measure the neutrino mass hierarchy. 

4.2 Measurement of the neutrino mass hierarchy 
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Figure 6: The gain fractions G (left panel) and G (right panel) for the mass hierarchy 
(MH) defined as in Sect. 2, i.e. the probability that the new facility will identify a given 
hierarchy at 5a provided that T2K+N0i'A have seen or not a hint of it at 90%CL irre- 
spectively of the value of 6. In both cases results are shown as a function of the year when 
T2K+N0vA stop taking data. 
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Fig. 6 presents the gain fractions G and G for the neutrino mass hierarchy 
discovery potential at So" under the hypothesis of a positive (left panel) or negative 
(right panel) result at T2K+N0t/A at 3ct, as a function of the T2K+N0z/A total 
running time. Consider first the gain fractions for NFIO and LBNE-33 kton. These 
two facilities have a very long baseline and sufficient statistics to measure the true 
neutrino mass hierarchy at Scr independently from the outcome of T2K+N0i/A. We 
indeed see that G — G —1 lor both of them. 

On the other hand, for phase I of the LBNE experiment, LBNE-10 kton, we have 
G ~ 1 whereas G ~ 0.85 depending on the outcome of T2K+N0z/A. The reason for 
this is that both T2K+N0z/A and LBNE-10 kton are more sensitive to the mass 
hierarchy for the same values of 5. Therefore, if T2K+N0z/A are not able to obtain 
a hint for the true mass hierarchy at 3a, then it is likely that 5 lies in the small region 
of the parameter space where LBNE-10 kton cannot reach the ha level either. 

Finally, the only facility whose gain fractions are strongly affected by the results 
of T2K-|-N0i/A is T2HK, whose very short basehne makes it very difficult to measure 
the hierarchy by itself. We see in Fig. 6 (left) that G increases with T2K-|-N0z/A 
running time, reaching a plateau around G ~ 0.6 by 2017. Conversely, if T2K-|-N0i/A 
are not able to find a 3(7 CL signal for the mass hierarchy with increasing exposure 
then the addition of T2HK will only allow a Scr discovery with a ~ 5% probability. 

5. Conclusions 

In this paper we have introduced a new method to compare the physics reach of 
future neutrino oscillation facilities, taking into account the impact of the outcome 
of present long baseline experiments (T2K and NOz/A). We have compared the new 
facilities on the basis of their capability to discover CP violation in the leptonic sector 
and to pin down the correct neutrino mass hierarchy, two of the missing inputs to 
complete our phenomenological understanding of lepton mixing. 

In order to compare facilities, we have first computed the probability P{B\5) that 
T2K-|-N0i/A find a 90% CL hint for CP violation or a Scr indication for the mass 
hierarchy for a given 5. We have done this by simulating 1000 possible outcomes 
of T2K-|-N0i/A by randomly generating realizations of the experiment properly dis- 
tributed around the theoretical expected values for fixed leptonic mixing parameters 
(including 5 and the neutrino mass hierarchy). We have then computed, apply- 
ing the same method, the probability P{A\5) with which one of the new facilities, 
in combination with T2K+N0z/A, is able to discover CP violation or the neutrino 
mass hierarchy. Eventually, we have computed the joint probability P{A,B\6) that 
T2K-|-N0z/A provides a hint for CP violation (or measure the mass hierarchy) at a 
given CL and that the combined data from T2K-|-N0i/A and the new facility make a 
discovery. The joint probability P{A, B\S) is a way to measure the complementarity 
of the new facility with T2K-|-N0i/A. Indeed, for the values of 5 where this proba- 
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bility is large, there is a significant probability that CP violation will be missed by 
T2K+N0z/A but still be discovered at the new facility. 

After having evaluated the degree of complementarity between T2K+N0z^A and 
the new facilities, we have introduced a very condensed performance estimator, the 
gain fraction G. This is the conditional probability P{A\B) = P{A, B)/P{B) that 
the new facility discovers CP violation (or the neutrino mass hierarchy) at 5a given 
that T2K+N0i/A have been able (G) or unable (G) to see a hint of it at the 90% 
CL (3(7), regardless of the value of 5. While the joint probability quantifies the 
complementarity between a new facility and T2K+N0z/A, the conditional probabilty 
allows to make well-defined statistical statements on how the performance of a new 
facility will be affected by the different possible outcomes of T2K+N0z/A. It is par- 
ticularly interesting to see how G and G for a given new facility scale with increasing 
T2K-|-N0i/A exposure, i.e. with increasing T2K-|-N0i/A running time. This is of 
extreme interest in a moment in which a decision on which facility should be built 
must be taken, since it helps answering the questions: What would happen if we keep 
running T2K-|-N0i/A beyond their scheduled data taking period? Is that helpful to 
improve the chances of new facilities to achieve their goals? Can the results from 
T2K-|-N0i/A help to decide what is the optimal next step? 

Armed with these concepts, we have compared three new facilities: T2HK, LBNE 
(in two versions, with a 10 kton and a 33 kton mass detector) and the 10 GeV 
Neutrino Factory (NFIO). The outcome of our analysis shows that: 

1. NFIO is, as expected, the facility which shows the most improvement over 
T2K-FN0i/A. For example, the joint probability P{A, B\5) {i.e., the probability 
that T2K-FN0i/A do not find a 90% CL hint for CP violation and that the 
combination with NFIO provides a 5a discovery) shows high peaks located at 
S ~ ±30°, ±150°, corresponding to regions of the S parameter space that can be 
tested by NFIO but not by T2K±N0z/A. On the other hand, peaks for T2HK 
(LBNE-33 kton) only appear for 6 = —90° (5 = 90°), i.e. at the same values 
for which T2K+N0t/A have the highest chances to discover CP violation. If 
we keep running T2K+N0i/A until 2030, the joint probabihties for these two 
facilities vanish, showing that the areas of the parameter space explored by 
T2HK and LBNE are less complementary to those probed by T2K+N0i/A. 

2. In case of a positive hint for CP violation at T2K±N0z/A at the 90% CL, 
NFIO has a 90% chance to improve the statistical significance of the hint up 
to 5(7. The gain fraction reaches a plateau very early, showing that running 
T2K+N0z/A beyond, say, 2020 will not increase the gain fraction of NFIO. 
The super-beam with highest chance to discover CP violation at 5a in case of 
a positive hint at T2K+N0i/A is T2HK. In this case, G also reaches a plateau 
for a rather low T2K+N0i/A exposure, and running T2K+N0i/A beyond 2020 
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does not increase their combined discovery potential significantly. The perfor- 
mance of T2HK, however, is strongly affected by our present ignorance of the 
neutrino mass hierarchy. If the hierarchy is still unknown at the time of T2HK 
data taking, then G ~ 0.45. On the other hand, if some other experiment is 
able to measure the hierarchy, we see that G increases significantly to ~ 0.75. 
Eventually, LBNE is the only facility whose gain fraction increases linearly with 
T2K-|-N0i/A exposure. In this case, keeping both experiments running adds 
information to the new facility data, reaching G ~ 0.45 if T2K+N0i/A keep 
running until 2030 (G ~ 0.1 when the 10 kton version is considered). 

3. In case T2K+N0i/A are not able to detect CP violation at the 90% CL, all 
of the considered super-beams will have a very low discovery probability, even 
though for low T2K-|-N0i/A exposure (say, until 2020) we can still see some 
hierarchy in the performances of T2HK and LBNE-33 kton (for which we get 
G ~ 0.2 and G ~ 0.1, respectively). If we keep running T2K+N0z/A until 2030 
with still no hint of CP violation, then the chances that one of these facilities 
will be able to discover CP violation at Sa will all be below 10%. On the other 
hand, the gain fraction G at NFIO ranges from G ~ 0.6 in 2020 to G ~ 0.35 
in 2030. This is the only facility with a non-negligible probability to discover 
CP violation in the leptonic sector if no hint of it is present even after running 
T2K+N0i/A until 2030. 

4. The discovery of the true neutrino mass hierarchy at 5o" is not a good observable 
to compare the performances of new facilities, as all of them (with the exception 
of T2HK) have been designed to maximize their probability to achieve this goal. 
However, it should be stated that if T2K-|-N0z/A detect a 90% CL hint on the 
true hierarchy by 2020, then the chances that T2HK can confirm such hint at 
5(7 are as high as 60% (all other facihties have G = 1 in this case). On the 
other hand, if T2K-|-N0i/A have not seen a preferred hierarchy by 2020, then 
the chances that T2HK can detect the true hierarchy at 5cr fall to less than 
5%. In this case, while NFIO and LBNE-33 kton still show (5 = 1, the smaller 
version of LBNE will have a very good chance to discover the neutrino mass 
hierarchy but it will not be guaranteed (with G ranging from 0.85 to 0.90 if 
T2K-|-N0i/A are kept running until 2020 or 2030, respectively). 
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